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Abstract 



O 

\ In this paper within the scope of FRW cosmology for A; = 0, ±1 we study the dynamics of the 



universe for a cosmological model with a scalar field nonminimally coupled to matter. By best- 



. fitting the model parameters with the observational data for the direct interaction between the 

! dark sectors in the model we obtain observational constraints on cosmological parameters. The 

■ result shows that with the best fitted model parameters, only in flat universe, the phantom crossing 

occurs twice in the past and once in the future, whereas no crossing occurs for open and closed 
models of the universe. 
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I. INTRODUCTION 



Recently, the observations of high redshift type la supernovae and the surveys of clusters 
of galaxies jl|-j4| reveal the universe accelerating expansion and that the density of matter 
is very much less than the critical density. Also the observations of Cosmic Microwave 
Background (CMB) anisotropies indicate that the universe is fiat and the total energy density 
is very close to the critical one jsl- 

The above observational data properly complete each other and point out that the dark 
energy (DE) is the dominant component of the present universe which occupies about %73 
of the energy of our universe, while dark matter (DM) occupies %23, and the usual baryonic 
matter about %4. There are prominent candidates for DE such as the cosmological constant 



a. 



a dynamically evolving scalar field ( like quintessence) js, 9| or phantom (field with 



negative energy) lOj that explain the cosmic accelerating expansion. Meanwhi 
celerating expansion of universe can also be described through modified gravity 
cosmology and so on 



e, the ac- 
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brane 



18|. 



To explain the early and late time acceleration of the universe, it is most often the case 
that such fields interact with matter; directly due to a matter Lagrangian coupling, indire ctly 
through a coupling to the Ricci scalar or as the result of quantum loop corrections 19|- 24 1. 
If the scalar field self-interactions are negligible, then the experimental bounds on such a field 
are very strong; requiring it to either couple to matter much more weakly than gravity does, 
or to be very heavy 25]-j28|. In some models, coupling of the scalar field to matter should 
be tuned to extremely too small values in order not to be conflict with the Equivalence 
Principal {29]. In the proposed models, the scalar field couples to matter with gravitational 



strength, in harmony with general expectations from string theory [30|]-32|. Then, the 
scalar field is permitted to couple to matter much more strongly than gravity does, and 
yet still satisfies the current experimental and observational constraints. The cosmological 
value of such a field evolves over Hubble time-scales and could potentially cause the late-time 
acceleration of our Universe 33[. The crucial feature that these models possess are that the 



mass of the scalar field depends on the local background matter density. Whi 



e the idea of 



a density- dependent mass term is not new [3^-[39|, the work presented in [32|, |33| is novel 
in that the scalar field can couple directly to matter with gravitational strength. Since the 
scalar field mimics the background (radiation/matter) matter field, subdominant for most 
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of the evolution history except at late times when it becomes dominant, it may be regarded 



cosmological tracker field 



40|. 



In this paper, we study the detailed evolution of a cosmological model that an scalar field 
nonminimally coupled to the matter lagrangian. We assume FRW cosmological model with 
open, closed and open background geometries. In Sec. 2, we present the model in terms of 
dimensionless dynamical variables. In Sec. 3, by best-fitting the model with the most recent 
observational data, the Type la supernovea (SNe Ia)for distance modulus and applying 
method, new constraints on the model parameters are revealed. Sec. 4 devoted to the 
cosmological test of the model by examining the behavior of the effective EoS parameter 
and also reconstructing the model evolution from the solution. The interaction model is 
presented in Sec. 5. We summarize our results in section 6. 



II. THE MODEL 



In this section we consider an action of gravity in the presence of a scalar field nonmini- 
mally coupled to the matter lagrangian, 

S = J - ^<7'^^V^</.V.</. + 1^(0) + f{<jy)L^]V^dx\ (1) 

where R is Ricci scalar, G is the newtonian constant gravity and (f>{x^) is the scalar field 
with a potential V{(f)). Unlike the usual Einstein- Hilbert action, the matter Lagrangian 
is modified as f{(p)Lm, where f{(p) is an analytic function of the scalar field. This last term 
in Lagrangian brings about the nonminimal interaction between the matter and scalar field. 
The variation of action ([T]) with respect to the metric tensor components in FRW cosmology 
yields the field equations, 

3(H' + ^) = p^f + U' + V{<P), (2) 

2 

2H + 3H' + ^ = -U' + V{<P), (3) 

where we put 8irG = c = h = 1 and assume a perfect fiuid with pm = ^Pm- In the above 
equations, for fiat, closed and open models of the universe. A; = 0, ±1, respectively. We also 
assumed that he energy density pm stands for the contribution from the cold dark matter 
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to the energy density, 7 = 0. Also variation of the action ([T]) with respect to the scalar field 
(f) provides the wave equation for the scalar field as 

4> + 3H^=-V' -Prnf, (4) 

where prime indicated differentiation with respect to (p. One can easily find that both 
components -the scalar field (p and the coupling with CDM, Pmf- do not conserve 
separately but that they interact through a term Q according to 

(pj) + sHpmf = g, (5) 



p^ + 3H{1 + 



(6) 



where Q = Pmf'4' is the energy exchange term derived from the model. We also define the 



ratio between the energy densities as r = 

In the following we present the structure of the dynamical system by introducing the 
dimensionless variables: 

^2 



g^2' -""Z 3^2' 3^2' — 3^2^2' 

where "mf" means coupling of the scalar function / with matter field. 

We consider that both scalar function and scalar potential behave exponentially as 
= /oexp(a0) and V{(/)) = Voexp(/30) where a and (3 are dimensionless constants 
characterizing the slope of potential V{(j)) and The cosmological models with such 

exponential functions have been known lead to interesting physics in a variety of context, 
ranging from existence of accelerated expansions 4l| to cosmological scaling solutions 42 1- 



Pmf 



V 



k 



(7) 



431]. In particular, the exponential forms of and V{(j)) are motivated by a variety of 



cosmological models 



44j and also from stability considerations 45|]. By using equations 



the evolution equations of these variables become, 
dill 



dN 

dQmf 

dN 

dfly 

dQk 

In 
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H 



mf 



mf 



mfi^ 



-2r2fc — 2r2fc 



H 



(8) 

(9) 
(10) 

(11) 



where N = ln{a) and 

= _l[3n^^ + 6n^]+fifc, (12) 

^ = Venj + + a^imf. (13) 

In term of the new dynamical variables, also the Friedmann constraint is 

nk + i = n^f + nv + n^. (14) 

One can write the effective deceleration and EoS parameters for the model and also the EoS 
parameter for the scalar field in terms of new variables as 

q = -l-Qk + SQ^ + ^Qmf, (15) 

By imposing the constraint f|T^ the equations fl8l)- f|TTl) reduce to three equations for the 
dynamical variables ^^^(0), f2m/(0) and f2y(0) . 

In the next section we solve the above equations by best fitting the model parameters 
a, P and initial conditions ^^^(0), r2mf (0), fiy(O) and H{0) with the observational data for 
distance modulus using the method. The advantage of simultaneously solving the system 
of equations and best fitting the model parameters is that the solutions become physically 
meaningful and observationally favored. 



III. COSMOLOGICAL CONSTRAINTS 



The difference between the absolute and apparent luminosity of a distance object is given 
by, fj.{z) = 25 + 5 logj^o '^^(-2^) where the luminosity distance quantity, d^i^z) is given by 

'^-w-'i-^^^f-ii)- '''' 

With numerical computation, we solve the system of dynamical equations for Q^, Q^, Q^f 
and Qy While best fitting the model parameters and initial conditions with the most recent 
observational data, the Type la supernovea (SNe la), in order to accomplish the mission, 
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we need the following two auxiliary equations for the luminosity distance and the hubble 
parameter 

^-<-^^)- 

To best fit the model for the parameters a, /3 and the initial conditions Q^{0), 17^/(0), 
Qv{^) and H{0) with the observational data, SNe la, we employe the method. From now 
on, we also assume that the matter field in the model is cold dark matter (CDM), 7 = 0. 
We constrain the parameters including the initial conditions by minimizing the function 
given as 

XsNe{(^^ /3; ^<^, ^m/, ^V, H\o) = 2 ' (^l) 

i=l « 

where the sum is over the SNe la data. In relation fl2Tl) . fif^'^ and are the distance 
modulus parameters obtained from our model and observation, respectively, and a is the 
estimated error of the fi"^'^. Table I shows the best fitted model parameters and initial 
conditions in all three cases, k = 0, ±1. 



TABLE I: Best-fitted model parameters 
model a (3 Q^{0) n^f{0) Qv{0) H{0) xLn 

k = 1.8 2.2 -0.4 -0.2 1.6 0.91 548.0229449 
k = 1 0.94 -0.07 0.7 -0.2 2.05 0.916 557.5901281 
k = -l 6.5 3.7 0.4 0.5 2.2 0.908 546.03629 



Figs. 1-3 show the constraints on the parameters a, (3 with the initial conditions on 
n^(0), rimfiO), ftv{0) and H{0) at the 68.3%, 95.4% and 99.7% confidence levels. The best 
fitted 2 dimensional likelihood the model parameters are also shown in figures 1-3. 
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FIG. 1: The best-fitted two dimension likelihood and confidence level for a, /3 in k = case 




FIG. 2: The best-fitted two dimension likelihood and confidence level for a, /3 in A; = +1 case 




FIG. 3: The best-fitted two dimension likelihood and confidence level for a, /3 in A; = — 1 case 
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FIG. 4: The best-fitted distance modulus fi{z) plotted as function of redshift in A; = 0, ±1 cases 

In Fig. 4, the distance modulus, fj.{z), in our model is fitted with the observational data 
for the parameters a, (3 and initial conditions ^2^(0), r2m/(0), fiv'(O) and i^(0), using 
method in all three cases k = 0, ±1. 

IV. RESULTS AND DISCUSSION 

In order to realize the behavior of the universe and its dynamics we need to verify the 
best-fitted model parameters. We have already best fitted our model with the current 
observational data for the distance modulus. The cosmological parameters analytically 
and/or numerically have been investigated by many authors for variety of cosmological 
models. The effective EoS and deceleration parameters, plus scalar field EoS parameter are 
given by the equations f|T5l) - f|T7j) . 

In Fig. 6, these parameters are shown for all three models of the universe. A; = 0, ±1. All 
the trajectories shown in these graphs are best fitted for the model parameters a and (3 and 
also initial conditions. From Fig. 6) left for scalar field EoS parameter one observes that in 
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FIG. 5: The best-fitted effective —^mf and, plotted as function of redshift in /c = 0, ±1 cases 

none of the closed, open and flat models of the universe the phantom crossing occur in the 
past and future. In k = case. The EoS parameter for the scalar field, u^, become tangent 
to the phantom divide line in the near past. For effective EoS parameter. Fig. 6) middle, 
phantom crossing occurs twice in the past and once in the future in flat model of the universe. 
The phantom crossing can also can be verified from Fig. 5, for A; = three times whenever 
= —Qmf for / < and for A; = ±1 never occur, as the conditions are not satisfied. The 
graph also shows in each case of open and closed universe, the behavior of EoS parameters 
are very similar to each other. The graphs for both EoS parameters in case of open universe 
also show a singularity in the near past. The deceleration parameter for all three cases 
exhibits a universe transition from deceleration state to acceleration one in the near past. 
In case of open universe the graph for the deceleration parameter approaches minus infinity 
in the near future, a phantom superacceleration that leads to Big Rip singularity. 

The dynamics of the interacting term, Q, and the ratio r are plotted in Fig. 7, for 
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FIG. 6: The best-fitted effective EoS parameter, scalar field EoS parameter and deceleration pa- 
rameter plotted as function of redshift in fc = 0, ±1 cases 

k = 0, ±1 cases. The sign of Q for different cosmological models depends on the sign of 
f{(p) and (p (Fig. 8). According to the plots in Fig.8, one finds that Q ^ for k = ±1, 
since (p is positive in both cases and / ^ respectively. It means that for k = —1 energy 
transfers from coupling pmj field to the scalar field fluid and visa versa for A; = +1. In flat 
universe, sign of the interaction term Q changes from negative in the past at 2; ~ 0.25 to 
positive in the near past and future, as f{(f)) is always negative and changes its sign at 
the same redshift. One may consider that the coupling pmf and scalar field act like dark 
energy and dark matter respectively. The ratio between densities for flat and closed model 
of the universe are negative whereas for the open universe is positive (Fig. 7). This is due 
to the presence of the function /((/>) in the relation that as can be seen in Fig. 8), for open 
universe is positive and for flat and closed universe is always negative. 

The evolutionary form of the best fitted coupling pmf, the reconstructed derivative of 
the scalar field 4> and potential V{(/)) are plotted in Fig. 8, for different cosmological models. 
The reconstructed scalar potential for open universe in Fig. 8 while vanishes in far past, 
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FIG. 7: The best-fitted interaction term Q and densities ratio r plotted as function of redshift in 
/c = 0, ±1 cases 
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FIG. 8: The best-fitted of Pm/, ^(0) and </> plotted as function of redshift in = 0, ±1 cases 



become negative in the late time and tends to minus infinity in the near future. In case of 
open and flat universes, tlie potential goes to -|-oo in the near future. 



In general, the interaction term Q can be an arbitrary function of cosmological parameters 
like hubble parameter and energy densities. In here, we assume three forms of Q as 

Q = amfHp„J, Q = o-4,Hp^, Q = aeffHpeff, (22) 
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FIG. 9: The best-fitted amf, cr^ and aeff parameters plotted as function of redshift in A; = 0, ±1 
cases 



where the interaction term is proportional to the three kinds of energy densities introduced 
in our model. From the best fitted interaction term derived in the model and with the 
observational data, one can describe the dynamics of the parameters, 0"^/? o"(/> and (Jeff- In 
Fig. 9, the parameters are plotted against redshift z for the cosmo logical scenarios k = 0, ±1. 
One common feature for all these parameters is that they are flattening in high redshifts. 
In closed universe scenario, the parameters vanish in future. The parameters are fiat in the 
past and only begin evolving in the near past. 



V. CONCLUSION 



In this paper, we have studied the evolution of the gravitational and scalar field in cos- 
mology in terms of dimensionless dynamical variables, in which a scalar field is nonminimally 
coupled to the matter lagrangian. The field equations are solved simultaneously by means 
of best fitting the model parameters with the data obtained from the observations of type la 
supernovae and utilizing the method for the distance modulus. This approach imposes 
observational constraints on the model parameters and yields physically more accurate so- 
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lutions and observationally more reliable results. We present the model for the cases of flat, 
closed and open types of the universe. 

The best fitted model parameters in flat universe for the effective EoS parameter show 
twice phantom crossing in the past and once in the future, whereas no crossing occurs for 
open and closed models of the universe. There are two necessary conditions need to be 
satisfied for phantom crossing, i.e. 1) f{(f>) < 0, 2) = —Qmf, that only in flat universe 
can be fulfilled. The graph for evolution of the effective EoS parameter in case of flat universe 
reveals a multiple transition from Ueff > — 1 (quintessence era) in the past to Wg// < — 1 
(phantom era) in the recent past, again to cjg// > — 1 very near past and finally to phantom 
era in future. The phantom crossing for A; = occurs within the range of the observationally 
accepted redshift z < 0.45. The best fitted deceleration parameter shows that in case of 
open universe the model predicts big rip in near future. 
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